arXiv: 1502.00347v2 [cond-mat.mes-hall] 15 May 2015 


Thermal vector potential theory of transport induced by temperature gradient 

Gen Tatara 

RIKEN Center for Emergent Matter Science (CEMS) 

2-1 Hirosawa, Wako, Saitama 351-0198, Japan 
(Dated: May 18, 2015) 

A microscopic formalism to calculate thermal transport coefficients is presented based on a thermal 
vector potential, whose time-derivative is related to a thermal force. The formalism is free from 
unphysical divergences reported to arise when Luttinger’s formalism is applied naively, because 
the equilibrium (‘diamagnetic’) currents are treated consistently. The mathematical structure for 
thermal transport coefficients are shown to be identical with the electric ones if the electric charge 
is replaced by energy. The results indicates that the thermal vector potential couples to energy 
current via the minimal coupling. 


Conversion of heat into electric and other currents and 
vice versa is of essential importance from the viewpoint 
of realizing devices with low energy consumption. Of 
recent particular interest is heat-induced spin transport 
in the field of spintronics, where spin current is expected 
to lead to novel mechanisms for information technology, 
and to devices with low-energy consumption due to the 
absence or weak Joule heating. 

A hot issue in spintronics is to use magnetic insulators, 
which are suitable for fast magnetization switching and 
low-loss signal transmission. Insulators have, however, a 
clear disadvantage that electric current cannot be used 
for its manipulation. Instead, temperature gradients be¬ 
come the most important driving force in inducing spin 
transport. To study thermally-induced spin transport 
theoretically, a microscopic formulation is necessary for 
full understanding and for quantitative predictions. A 
microscopic description is, however, not straightforward; 
temperature gradients and thermal forces are macro¬ 
scopic quantities arising after statistical averaging, and 
thus it is not obvious how to represent those effects in a 
microscopic quantum mechanical Hamiltonian. 

In 1964, Luttinger proposed a solution [lj. To describe 
the effect of temperature gradient, he introduced a scalar 
potential T, which he called a ‘gravitational’ potential, 
which couples to energy density of the system, £, via an 
interaction Hamiltonian, 

if L = J d 3 r'S>£. (1) 

Although the microscopic origin of the potential has not 
been addressed, he argued that to satisfy the Einstein re¬ 
lation the potential adjusts itself to balance the thermal 
force, resulting in an identity VT = in the thermal 
equilibrium. Owing to this trick, thermal transport coef¬ 
ficients can be calculated by linear response theory with 
respect to the field 4' without introducing the tempera¬ 
ture gradient in a microscopic Hamiltonian. Another ap¬ 
proach, based on the Landauer-Biittiker formalism, was 
presented by Butcher p]. 

Luttinger’s method has been applied to study vari¬ 
ous thermally-induced electron transports M , magnon 


transport 0 and thermally-induced torque (11]. It 
turned out, however, that naive application often leads 
to apparently wrong transport coefficients which diverge 
as T —> 0 0- 1 11 the case of the thermal Hall effect, 

the divergence was identified to be due to a wrong treat¬ 
ment of the equilibrium diamagnetic current induced by 
the applied magnetic field, and it was found that the 
physical Hall coefficient is obtained if one subtracts the 
equilibrium magnetization current before applying linear 
response theory 0. A similar problem was reported re¬ 
cently for thermally-driven spin-transfer torques [ul- 

In the case of electrically-driven transport, elimination 
of unphysical equilibrium contribution from transport co¬ 
efficients is guaranteed by U(l) gauge invariance, which 
represents charge conservation. In the presence of an 
electromagnetic vector potential, A, the physical electric 
current has two components, a paramagnetic current (the 
first term) and a diamagnetic current (the second term), 
as 3 = m (p) ““ f^n e A, where (p) is quantum average 
of the momentum density operator, n e is electron den¬ 
sity and e and m are electron’s charge and mass. The 
paramagnetic current contains an equilibrium contribu¬ 
tion arising from all the electrons below the Fermi level, 
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which turns out to be — n e A. This equilibrium contri¬ 
bution thus cancels perfectly with the diamagnetic con¬ 
tribution, leaving only the contribution from excitations 
in the transport coefficients 0 ■ Obviously, a consistent 
treatment of the two contributions is necessary for the 
cancellation of equilibrium contribution and for gauge- 
invariant physical results. If one uses, instead of a vector 
potential, a scalar potential to describe an conservative 
electric field, the role of the diamagnetic current is not 
clearly seen, and wrong results easily arise if an inconsis¬ 
tent treatment is employed. 

From those experiences in electrically-induced trans¬ 
port, the divergence in the thermally-induced transport 
described by Luttinger’s ’L is expected to be due to an 
incorrect treatment of the ‘diamagnetic’ contribution. [0 
What Qin et al. 0 showed is that the ‘diamagnetic’ effect 
is consistently taken account of if divergenceless magneti¬ 
zation current is included. If one could construct a vector 
potential representation of thermal effects, the ‘diamag- 
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netic’ effect would be treated consistently and straight¬ 
forwardly, since the ‘diamagnetic’ current is defined by 
the Hamiltonian and the vector potential without ambi¬ 
guity. 

Temperature gradients exert a statistical force propor¬ 
tional to VT, which is conservative, i.e., has no rota¬ 
tion component. Still, one may introduce a rotation¬ 
less vector potential to describe the force. In the case 
of classical charged particles described by a Hamiltonian 
H = < ' p ~ 2 ^' > + V( r ) (V is a scalar potential), the total 
force is F = — — W if V x A = 0. Any rota¬ 
tionless force can thus be represented by use of a vec¬ 
tor potential without introducing a scalar potential. In 
thermally-driven transport, a thermal force proportional 
to VT is represented by a vector potential which we call 
the thermal vector potential. 


The objective of this paper is to propose a formal¬ 
ism describing thermal effects by a thermal vector po¬ 
tential, and to demonstrate that the formalism works 
perfectly for a few simple cases of thermally-driven elec¬ 
tron and energy transport, without yielding unphysical 
divergences. Since the ‘charge’ to which the temperature 
gradient couples is energy, the thermal vector potential 
couples to the energy current density operator, jg. We 
first carry out a derivation of a thermal vector poten¬ 
tial form of the interaction Hamiltonian by looking for a 
Hamiltonian equivalent to the Luttinger’s Hamiltonian. 
Local thermal equilibrium is assumed. We then derive 
expressions for electric current and energy current by use 
of conservation laws, and identify the ‘diamagnetic’ cur¬ 
rents. It is shown that the ‘diamagnetic’ currents remove 
unphysical equilibrium contribution to transport coeffi¬ 
cients. The results satisfy the Wiedemann-Franz law. We 
shall demonstrate that the obtained expressions for the 
currents indicates the minimal coupling of the thermal 
vector potential. 


It was recently demonstrated by Shitade that a gauge 
theory of gravity constructed imposing local space-time 
translation symmetry contains a vector ‘gauge field’ as 
well as a scalar potential corresponding to that of Lut- 
tinger 14f. The model was applied to describe thermal 
transport of non-interacting electrons and Wiedemann- 
Franz law was shown to be satisfied. The origin of local 
translation symmetry in the context of thermal transport 
was not addressed to. 


In this paper, we use energy conservation law to derive 
a vector potential representation of thermal effects, with¬ 
out assuming invariance under local space-time transla¬ 
tion. We start with rewriting the Luttinger’s Hamilto¬ 
nian by use of continuity equation for operators £ and 
3 £, 


S = -V-j e . 


( 2 ) 


as 

H h (t) = f d 3 r f dt'j £ {t’) ■ Vf(r,l) (3) 

J J — OO 

where we used Gauss’s theorem assuming that no field 
exists at r —» oo. This expression is not of the form 
of an interaction between a vector potential and energy 
current because of the time-integration. We here look 
for a Hamiltonian Ha t which agrees with Eq. © when 
long time average is considered, namely, f_ dtHA T = 
fZ dtHjXt). The result is 


H At = - J d3r js ( r > t) ■ At (t) (4) 

where At(£) = df'VT^') [l5! is the thermal vector 
potential, which satisfies 


d t A T (r,t) = VT(r,f) = —. (5) 

The interaction Hamiltonian © is understood as rep¬ 
resenting the thermodynamic potential change when a 
static temperature gradient is applied. In fact, the rate 
of the change of the entropy (5) due to an energy current 
is 0 


/ i r vt 

d 3 r-V ■j e = -J d 3 rj s ■ ( 6 ) 

and this entropy change modifies the thermodynamic po¬ 
tential, E — TS — fiN (E is the internal energy and N 
is the electron number). The effective Hamiltonian de¬ 
scribing a DC thermal force, H$ = —TS : is therefore (to 
the linear order in VT) 

H s = ^J d 3 r f js{t')dt' ■ VT. (7) 


This reduces to Eq. © after the replacement VT —>• 

We now apply the thermal vector potential interaction, 
Eq. ©, to study thermal transport and demonstrate 
that the formalism works perfectly. We consider free 
electrons with a quadratic dispersion, described by the 
Hamiltonian Ho = fd 3 r£ 0 , where £q = 2^.(Vc' l ’)(Vc) — 
fie' c is the free electron energy density, is the chemical 
potential and and c are creation and annihilation op¬ 
erators of the electron, respectively. The energy current 
density is derived by use of the energy conservation law, 
Eq. ©. For free electrons, V ■ = — j-[H 0 , fo( T ')] I and 

the result is 1121 


,-(°) - 
J£ ~ 


ih 3 

(2m) 2 


(V 2 c t )Vc— (Vc f )(V 2 c) 


^■( 0 ) 

J 1 

e 


( 8 ) 


where j ^ -c t V c is the paramagnetic part of 

the electric current density. We focus on the uniform 
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component of the current considering the case of spatially 
uniform temperature gradient, which reads (V is system 
volume) 


( 0 ) 

J £,i ~ m \' 


V 




( 9 ) 


where T = i^A T j £ fcw 7*f f{ u )<t>ku agrees with the 
diamagnetic current (Eq. (1101) 1. The equilibrium (dia¬ 
magnetic) contribution is therefore eliminated from the 
physical thermally-induced electric current, obtaining 
ji = cr T E T ,i , where 


where fc is a wave vector and e k = — A* is the energy 

measured from the Fermi energy. 

We now apply this interaction to study thermally- 
driven longitudinal electron transport on a basis of di¬ 
agrammatic (Green’s function) formalism. Besides the 
interaction Hamiltonian, we need to take account of the 
‘diamagnetic’ current contribution proportional to At- 
We derive it by use of the charge conservation law, 
p + V • j = 0 (p is electric chage density) taking ac¬ 
count of thermal vector potential. Namely, we calculate 
a commutator, — j[[Ha T ’ c ^c] = V ■ j Ar , and derive the 
expression for j Ar . The result of the uniform component 
is 


Eji = — 


dA T (r,t ) 

at 


VT 

= - —, ( 15 ) 


Ji 


At = ~—A t * 


*-T,j 


where 




h 2 


lk c k c k 


7fc = e k$ij + 


( 10 ) 


( 11 ) 


is the thermal field, and gt 

6m 2 Vr 2 J2ku fcW'MIsjJ 4 (assuming rotational 
symmetry for fc). The low temperature behavior is seen 
by a series expansion (<&*,(<*;) = jc^l 4 ) 


/ dw/'(w)$*(w) =-* fc (0) - — (fc B T) 2 $"(0), (16) 

J -oo 6 

where 0(T 4 ) is neglected. Since e k = 0 on the Fermi 
surface, we have £ fc fc 2 efc<l>fc(0) = 0 (to the leading order 
of ^L_). We therefore see that or = 0(T 2 ) at T —> 0 
and the thermally-induced current vanishes at T = 0. 

We now study the thermally-driven Hall effect and 
show that the vector potential formulation does not lead 
to unphysical results like the one reported in the Lut- 
tinger’s scheme Q- We introduce the interaction with 
an electromagnetic vector potential, Ha = — f d 3 rA ■ j, 
to describe the effect of the applied magnetic field, where 


j = J (0) 


j Ar , j A = — ^-AAc being the diamag- 


As expected from Eq. the diagrammatic calcula¬ 
tion is carried out by a straightforward replacement of 
charge e in the electric field driven case [12| by energy 
ek ■ Including the interaction with the vector potential to 
the linear order, the DC paramagnetic current is 


J i 


.(o) _ eS 


= _ V e 3 
mV k 

kco 


- 2A Tj /Mlm[(<4J 2 ] 
( 12 ) 


where ^kikje k , A T = <j) ku] = 


Z^f d “ 


0 phui 


27r 

+ 1] _1 is the Fermi distribution function (/3 = 
(fcT) -1 , kg being the Boltzmann constant). The re¬ 
tarded and advanced Green’s functions for free electron 
are denoted by . and 9%, u = (9 T k ,J*, 

where r is the elastic lifetime. By use of ^f-(g kLU ) 2 = 
dkjPku) and integration by parts with respect to fc, we 
rewrite the last contribution using 


j^k, — at , Yku) — g k ,u 9k,u> 

Im denotes the imaginary part and /(w) = 


netic current of the electromagnetic origin. The elec¬ 
tromagnetic vector potential is treated as static but has 
finite wave vector, since its role here is to represent a 
static magnetic field, B = V x A. The thermal vector 
potential has an infinitesimal angular frequency (fl) and 
is spatially uniform. The Hall current is calculated to the 
lowest order, i.e., linear in both A and At- It turns out 
that the leading contribution is linear both in the angular 
frequency f l and in the wave vector q. 

The contributions to the paramagnetic part of the cur¬ 
rent, j are shown diagrammatically in Fig. ®a)(b). 
(‘Diamagnetic’ currents shown in Fig. [Qc), vanish, since 
A and At carry only either finite angular frequency or 
finite wave vector in the present description.) The con¬ 
tribution of Fig. Oja) is 


.(a) 

Ji = 


e 2 h 4 

m 3 V 


(V m A j )A T , l J2k i k j 'y l k m $ ( E!, ( 17 ) 


kco 


E< 


(ffL) 2 - (sD 5 


= -'52'Yk‘fiku 


to obtain 


,-( 0 ) 


Ji = 


efc v > 
mV ^ 

kco 


\ ij f (k-0 / / \2 -Aj 

—-—((frkcu) -Ji 


( 13 ) 


( 14 ) 


where $£> = Im [/»(<£, J 2 «£ >(j , + fi/MG^J 4 ]. 

Because of ‘diamagnetic’ current due to the thermal 
vector potential, j Ar , we have an interaction vertex, 
— fd 3 r A ■ j AT , linear in both A and At- The contri¬ 
bution shown in Fig. [Ub) arises from this interaction 
vertex. It is 


jf } 


: ^y{VmAj)A T ,l Y, hkmik® 

kco 


2 


( 18 ) 
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M 



phenomena leads straightforwardly to physical transport 
coefficients. The uniform contributions to the electric 
and energy current densities we have derived are 





FIG. 1. Diagramatic representation of the contributions to 
the thermal Hall effect. Solid, wavy and dotted lines denote 
the electron, thermal vector potential At and the electro¬ 
magnetic vector potential A, respectively. Diagrams (a) and 
(b) correspond to contributions of Eqs. (1171) and (1181) . respec¬ 
tively, while contributions of diagram (c) vanish. 


The total Hall current, j® = j^ + j( h \ is finally ob¬ 
tained as 


J?° = Qh(Et x B), (19) 


cH 1 ^—\ 

Ji = 77 / v 

m 1 / a —' 


j£,i — 


m V 

h l 
m V 


k L 


ki - eAt - A T ,j 


CfeCfe 


E 


k L 


he k - e^Aj - 7r,k A T,3 


c k c k- ( 22 ) 


The key for direct access to physical results in the present 
formalism is the particular relation between the interac¬ 
tion vertex and ‘diamagnetic’ contributions to currents, 
such as 7 ^ = d ki (kje k ) and = ^:[kj(e k ) 2 ]. We 
finally show that these identities indicate that the to¬ 
tal Hamiltonian, H , including the electric and thermal 
vector potentials are of the minimal form (to the second 
order in thermal vector potential) [Tij] , 


where Q H = 7ku k 2 e k ^ k ]. We see that the Hall 
current vanishes at T = 0 (see Eq. m), as is physi¬ 
cally required. The vector potential formalism applied 
straightforwardly therefore leads to the correct result, in 
sharp contrast to Luttinger’s ‘gravitational’ potential for¬ 
malism. 

For consistency of the vector potential formalism, we 
need to confirm that ‘diagmagnetic’ current arises also for 
the energy current. This is not a trivial issue, since we 
cannot invoke a gauge invariance concerning the energy 
current, in contrast to the case of electric current. In our 
scheme, ‘diamagnetic’ contribution is explored again by 
looking into the energy conservation law. In fact, includ¬ 
ing the thermal vector potential interaction, Eq. 0, in 
the left-hand side of Eq. (0 , we see that the energy cur¬ 
rent acquires a ‘diamagnetic’ contribution linear in At- 
After a straightforward calculation, its uniform compo¬ 
nent is obtained as 

3s A = ~~ A T,jy7^jT k c l k c k , ( 20 ) 

k 

where 

( 2 /i 2 \ 

7t,/s = ^k f ^k^ij T kikj J . (21) 

We see here that the matrix 7 j! k for this energy cur¬ 
rent correction satisfies 77 ? k = ~^-[kj(k k ) 2 ]. and thus 
cancellation of the unphysical equilibrium contribution 
occurs, in the same manner as the electric currents dis¬ 
cussed above. The total energy current density induced 
by the thermal vector potential is jg = -kVT, where 

K = 2Ft (^) 2 Efc M fc i fc l( e fc) 2 / , M(^) 2 - The coeffi - 

2 2 

cient satisfies the Wiedemann-Franz law, k = — ks Ta b, 
where ctb is the Boltzmann conductivity. 

We have confirmed that thermal vector potential for¬ 
malism applied for various thermally-induced transport 


ft 2 

H =—^(k-eA-e k - eA A T ) 2 clc k . (23) 

k 

In fact, as is easily checked, -jj- = m(ki—eAi—'y k Atj) 
to the linear order of vector potentials (4j- denotes 
a functional derivative), and thus the formal definition 
of current, j = —jj~, agrees with Eq. (B^l) . As for 
the energy current, it is formally defined by jgg = 

+ s^~c) d ) ■ B y use of the Heisenberg equa¬ 
tion of motion for c k and c' k , we see that this formal 
definition agrees with Eq. (l22l) . 

In the electromagnetic case, the minimal form is im¬ 
posed by a U(l) gauge invariance. For the thermal vec¬ 
tor potential, in contrast, there is no gauge invariance 
in the strict sense since the energy conservation arises 
from a global translational invariance with respect to 
time. [l9j Still, our analysis indicates that the minimal 
form emerges. This fact might be understood as due to 
a ‘gauge invariance’ as a result of the energy conserva¬ 
tion law. In fact, we have shown that the Luttinger’s T 
and the present At have the identical effect concerning 
steady state properties. In other words, we may assign 
a part of thermal force to T and the rest to At, so that 
yH- = V\H + At- Thus we have a ‘gauge invariance’ un¬ 
der a transformation 41 —> T — x and At a- At + Vy (x 
is a scalar function). Such a gauge transformation is gen¬ 
erally defined for a vector field coupling to a conserved 
current. 

The thermal vector potential formalism applies to ther¬ 
mal torque straightforwardly by replacing the electric 
charge in the derivation of Ref. [20( by energy. It is 
easy to check that the thermal torque vanishes at T = 0, 
as is physically required. 

To conclude, we have demonstrated that the calcula¬ 
tion of transport coefficients are straightforwardly carried 
out based on a vector potential formalism. The formalism 
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would apply straightforwardly to the case of multi-bands 
and to the interacting cases. 
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